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Целью работы является исследование скорости сходимости 
четырехслойной итерационной схемы. Рассматривается 
задача нахождения приближенного решения линейного 
операторного уравнения Аи = } Для решения такой задачи 
используются двухслойные и трехслойные итерационные 
методы. При этом трехслойные методы сопряженных 
направлений сходятся значительно быстрее, чем двухслойные 
градиентные методы. Задача исследования —Й установить, 
имеет ли четырехслойная схема преимущество в скорости 
сходимости по сравнению с трехслойной схемой. Для этого 
приводится четырехслойная итерационная схема решения 
сеточных уравнений, и рассчитываются ее параметры. 
Доказано, что  четырехслойная  итерационная схема 
вариационного типа для решения сеточных уравнений 
выражается к трехслойной схеме. 
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ТБе уогКк оБесвуе 1$ ю зш4у Ше юпг-1ауег зсВете сопуегеепсе 
гае. Тве ргоет оЁ Нп@ те ап арргохипае зоаНоп ю Ше Ппеаг 
орегафог едпайоп Аи = } 1$ сопз14еге4. Т\ууо-1ауег ап4 Шгее-1ауег 
Цегануе те@о4$ аге изе4 1ю зо1уе 1$ ргоМеш. А фаё, Фе @гее- 
1ауег сопасае АшесНоп$ тефо@$ сопуегое Ёазег Фап Фе 1%0- 
1ауег отаФеп тефоа$. ТЬБе гезеагсь ргоМет 15 ю езмаБИзВ 
урвешег фе гопг-1ауег зсВете Баз а зрее4 а4уап'асе аз сотрагеа 
(0 Ше Шгее-Йауег зсВете. ТБе оиг-ауег зсВете 1$ сопзгасеа, 
апа Из рагатейег$ аге са1сша(еа юг 15 ригробе. И 1$ ргоуе4 фае 
Фе опг-Йауег Иегануе зсБете оЁ а уапанопа! 1уре Юг зоушз 


Воие-а!егепсе едиаНопз 4о\/п$ 10 Ше @гее-1ауег зсБете. 


Кеууогд5: Нппе-Ч@егепсе едиаНоп$, гее-1ауег зспете, ючг- 
1Лауег зсБеште, уапайопа! тефоа$. 


Введение. Большинство прикладных задач таких, как задача транспорта веществ [1-3], гидродинамики 
мелководных водоемов [4—5], аэродинамики [6-7], динамики популяций [8] и других, сводятся к решению системы 


линейных алгебраических уравнений (СЛАУ). Для решения таких систем уравнений используются двух- и 


трехслойные итерационные схемы. 


Рассмотрим задачу нахождения приближенного решения линейного операторного уравнения [9]. 
Аи=Х, (1) 


где А — симметричный положительно определенный оператор, действующий в вещественном гильбертовом 


пространстве Н. 


Для увеличения скорости сходимости вместо двухслойных итерационных методов используются трехслойные 
итерационные методы. Эти методы исследованы в работе [10]. Ниже приведено исследование четырехслойной 
итерационной схемы. Условия устойчивости такой схемы получены в работе [11]. 

Четырехслойная итерационная схема решения сеточных уравнений имеет вид 


ВУ = Вьн (В- ин А) ук + (1- ау) Вр + (ол — Ве ) ВУр-> +Вьыткы (2) 
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для к = 2,3...., Ву, = В(В-т,А) у +, Ву› = В» (В-т›А) у, +(1-а.› ) Вуз +В-т>/, %ЕН. 
Необходимо найти параметры {тк} : [© к} и {Ви} ‚ при которых норма эквивалентной погрешности 
х; = у, -и была бы минимальной для любого (. 
Расчет параметров схемы. Перепишем (2) в виде 
В (Увы + (вы —1) ука + (Вы — бы } Ук Ве — У 


ТАН 





НЯ 
Действительно, для уравнения погрешности схемы (2) 
хр + (9—1) жк (Вы — и ) к 
Вен 
Или ху = Вик (Е С) хь + (1- а) м + (01 Ван) 2. (3) 


Для минимизации нормы х, вН (п>1) необходимо и достаточно, чтобы 





ж=-тиСхк, С= ВТА. 


[инь сх,]=0, =, @) 
ИЛИ (янь СХ, | = Ри (Схь,Сх,) =0. 
При К =] получим (хил, Схь ) = Вити (Сх›Схь)=0. 
Из (3), (4) следует 
(Сом ) = Вин (Св, (Е — аб) к ) + (1 — Чт (Сож ) + 
+ (+1 — Выы )(Схь-2, 2), 
(Сжльжн) = Вы (Сж-ь(Е пы) жь + (1-х )(Сжльж)+ 
+ (а — Вкл} (Схь-ьХь-2) › 
(С Же) = Ву (СХ › (Е С) ) с (1 О (С 1) + (0 —Вьнь )(Схь, №) . 
Запишем систему для расчета {тк} з [д } и {В} 
скан (СХ, Схь )+ (ал — Вкл) (Сжь-2, к ) = 0, 
—срыВьы (Си, Схь )+ (1 оу) (Сжль жи) = 0, 
(Ся, )-тьи (Сы Сх, )=0. 
Преобразуем систему уравнений 


(Ся, №2) 
Ск, Ху ) м (Со, Сх ) 





Вы = 1 ( 


(бо: Хо СЖ Сх, ) 





+ (Са ) тии (С: СХ, ) 1 (Се, ХА ) = (Схьь Ж), 
и (СХ 5х ) 
ы (СХ ‚С ) | 


Введем обозначение 


(Сож) 
(Схко, Хк2 ) Ан (С Схь ) 


Г (Сх ый) И {о т. № 
Тин = ‚ Чи = › 
(Сл СХ ) (Ся, 1 ) + фин (Ся, СХ ) 


Винт = О афкн. 


Преобразуем выражение (3) 
Схь-2 = (же + Вь-ихь-2 + (1-1) жь-з + (ви — В) 4 (сы } 


Запишем выражение (Ся, Сх, ) с учетом полученного выражения 


(Сх2»Сх,)= (5 Ве + (1-01) жк-з + (1 Вы) ма (ть Вк- С ) =0. 


фин = 





‚ тогда 
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Таким образом, получим ф,.; =1. Следовательно, Ву.1 = ох. 

В игоге выражение (3) преобразуется к виду 

Жен = 0 (Е ии С) жи + (-би )хьф. 

Выводы. В итоге получили, что х„., зависит только от х„, Х_| и не зависит от х„, п=0,К-2. Другими 
словами, доказано, что четырехслойная итерационная схема решения сеточных уравнений преобразуется к 


трехслойной схеме, поэтому использование первой не дает увеличения скорости сходимости по сравнению со второй. 
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